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We present a Schwinger-boson approach to the Heisenberg model with 
Dzyaloshinskii-Moriya interaction. We write the anisotropic interactions in 
terms of Schwinger bosons keeping the correct symmetries present in the spin 
representation, which allows us to perform a conserving mean-field approxi- 
mation. Unlike previous studies of this model by linear spin-wave theory, our 
approach takes into account magnon-magnon interactions and includes the 
effects of three-boson terms characteristic of noncollinear phases. The results 
reproduce the linear spin-wave predictions in the semiclassical large-5 limit, 
and show a small renormalization in the strong quantum limit S = 1/2. 

For the sake of definiteness, we specialize the calculations for the pat- 
tern of Moriya vectors corresponding to the orthorhombic phase in La2Cu04, 
and give a fairly detailed account of the behavior of ground-state energy, 
anisotropy gap, and net ferromagnetic moment. In the last part of this work 
we generalize our approach to describe the geometry of the intermediate phase 
in La2-xNda;Cu04, and discuss the effects of including nondegenerate 2p z oxy- 
gen orbitals in the calculations. 
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I. INTRODUCTION 



Dzyaloshinskii |1J was the first to point out that weak ferromagnetism in antiferromag- 
netic compounds can be explained by an antisymmetric spin-spin interaction. The mi- 
croscopic basis for this theory was later given by Moriya, @ who extended Anderson's 
superexchange theory |3[ to include the spin-orbit interactions. He considered that only the 
antisymmetric part of the anisotropic spin-spin interactions he derived was relevant because 
it is linear in the spin-orbit coupling, while the symmetric anisotropies are of second order. 
However, in a recent work Shekhtman et al. ||] showed that the symmetric anisotropic term 
cannot be neglected, and that the whole two-spin anisotropic interaction can be mapped 
onto the isotropic Heisenberg interaction by a bond- depending rotation of the original spins. 
This surprising result means that the anisotropy coming from the spin-orbit coupling does 
not lift by itself the ground-state degeneracy, unless the required local transformations to 
map the entire Hamiltonian to the isotropic one are not all consistent. Using the map- 
ping to the isotropic model these authors showed that infinitely many states with different 
ferromagnetic moments are degenerate with the purely antiferromagnetic one. Then, they 
concluded that some kind of frustration of the required local transformations is a necessary 
ingredient to have a definite ferromagnetic moment in the ground-state, a condition that 
had been previously overlooked. 

The results of Ref. Q are relevant to explain the small ferromagnetic moment observed 
in the Cu0 2 planes of undoped La 2 Cu04. This compound undergoes a structural phase 
transition from tetragonal to orthorhombic symmetry at T ~ 500 K. || In the orthorhombic 
phase the CuC>6 octahedra forming each CuO layer tilt in a staggered pattern about the 
(110) axis, and this distortion combined with the Dzyaloshinskii-Moriya (DM) interactions 
induces a weak ferromagnetic moment in each layer. Coffey et al. determined the 
DM vectors in both phases from symmetry analysis and detailed microscopic calculations. 
They emphasized that the existence of a net ferromagnetic moment requires an alternating 
pattern of Moriya vectors that change in sign from one bond to the next. However, only 
the consideration of the hidden symmetry discussed in [Q allows one to show that only 
for the orthorhombic symmetries the corresponding pattern leads to the observed small 
ferromagnetic moment, while those of the tetragonal phase do not allow it. These results 
where later generalized by Bonesteel M to describe the intermediate phase observed in 
La2- x Nd x Cu04, where the oxygen octahedra rotate around a general axis in the CuC>2 
plane. [[K|] 

In this work we focus on the study of the Heisenberg Hamiltonian with DM interac- 
tion in the strong quantum limit 5* = 1/2. The same Hamiltonian was previously studied 
by semiclassical linear spin- wave theory in Refs. [P|, |TT|] , with the emphasis put on magnon 
dispersion relations. However, since the classical ground state is not collinear, relevant 
three-boson terms coming from the Holstein-Primakoff representation of spins have been 
ignored. In this case we extend the Schwinger-boson theory discussed in [|12| to consider the 
anisotropic DM terms, and include magnon-magnon interactions in a saddle-point evalua- 
tion of the resulting bosonic Hamiltonian. We give a fairly detailed account of the energy, 
ferromagnetic moment, and gap behaviors as a function of the size and relative orientation 
of Moriya vectors, and compare these results with similar predictions from linear spin-wave 
theory. We will first specialize the calculations for the pattern of microscopic Moriya vectors 
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corresponding to the orthorhombic phase in La2Cu04. Then, in the last part of this work, 
we will generalize our approach to describe the intermediate phase in La 2 _ a; Nd :r Cu04. 

The layout of this work is as follows. In Sec. 2 we discuss the representation by Schwinger 
bosons of the anisotropic terms coming from the DM interaction. In Sec. 3 we describe 
the classical ground state of the model, and interpret its structure as a spiral order in a 
topologically-equivalent decorated lattice. This allows us to make contact with the theory 
developed in [T3]. In Sec. 4 we sketch the calculation of quantum fluctuations above this 
classical ground state. We give here explicit expressions in terms of the order parameters for 
the quasiparticle dispersion relation, ground-state energy, etc. In Sec. 5 we solve numerically 
the consistency equations and present the results for relevant quantities. In addition, in Sec. 
6 we generalize the calculations to describe the intermediate phase in L^-xNd^CuC^ and 
discuss the effects of including the oxygen 2p z orbitals. |TTl , |l"3"f Finally, we end up with some 
conclusions. 



II. BOSON REPRESENTATION OF THE DZYALOSHINSKII-MORIYA 

INTERACTION 

In its simplest derivation, the Heisenberg Hamiltonian including the DM interaction 
can be obtained starting from the one-electron Hamiltonian 

^i- ei = E e ^W+ E ci(te ieZ "-°U'Ci*' , (1) 

i,cr (ij},crcr' 

where i,j are nearest neighbor sites, and the term nondiagonal in spin space takes into 
account the spin-orbit interaction. We assume nondegenerate ground-state orbitals, and t 
and 9 are taken independent of the bond (ij) for simplicity. The unit vector = — dji. 
Then, by using second-order perturbation theory in the presence of an on-site repulsion U 
one obtains the interaction between spins at sites i,j as 

/•:,, = JSi.Sj + i^.iS; x s ; ) + s,,r, ; ,s ; . (2) 

In this expression the superexchange interaction J = (At 2 /U) cos 2 9. The last two terms 
of @ are the contributions coming from the spin-orbit interaction of the original electrons. 
They contain the antisymmetric Moriya vector ~Df- = 2Jrdij (r = tan#), and the symmetric 
anisotropy tensor T. In fact, (|J) can be most simply written as @ 

Eij = J (Si.Sj + 2rd ij -.(S i x Sj) + r 2 [2(^.3^^) - S ?; .S 

= J Sj.Sj , (3) 
where Jo — At 2 /U. The spins S^,S^ are obtained by rotating the original ones around the 

S' k = (1 - cos9 k )(d ij .S k )d ij + cos9 k S k - sin0 fe S k x d ij( - (k = i,j), (4) 
with the rotation angles given by 9i = 9 = —9j. 
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At this point we introduce the Schwinger representation S' = a't.er.a', where the bosonic 
spinor a't = (a+,an). Then, the rotational invariant structure in (|3]) can be written in the 
usual way as 

S'i-S'j =: 4'4 : —A^A'^ , (5) 

in terms of the SU(2) singlets 

B 'i0 = 9 J2 a t a 'ja > A 'ij = 9 aa 'ia a 'j,~v ■ 
Z a Z cr 

In (|5]) the colons denote a normal-ordered product of Bose operators. 
From standard rotation properties we have: 

S' = K 6 S = aJ.(R e a).a = (a + z4).«7.(Z^a) = a' f .a.a' , 

where IZg is the three-dimensional matrix that rotates an angle 9 around the axis, and 
Uo = exp[— i(6/2)dij.a) is the corresponding SU(2) matrix. In this way we have a' = Z^a, 
and this relation can be replaced in the definitions of the A', B' singlets in order to express 
them in terms of the unrotated spinors: 

B% = cos QB\ 3 + sin QC\ 3 , A\- = cos 6A tJ - sin 6D l3 . (6) 

Here, in addition to the singlets A, B we have defined 

^Ij - 2 ^ a \a{id-ij- a )™ ,a jv' j Aj = 2 a itT ((T v d ii .o) tT<T ia ja > . 

acr' a, a' 

By replacing (^),(^) in (|3|) we can identify the correct expressions in terms of Schwinger 
bosons of the spin interactions: 

Sj.Sj =: BjjBij : — A^A^j 
dy.(Si x S,) = i(: + : +4Ai + (7) 

2(dij.Si)(dij.Sj) — Sj.Sj =: CjjCij : —D^Dij . 

In this way the DM terms are given as products of bilinear boson operators that have the 
same rotational properties of the original spin interactions. This allows us to perform a 
conserving mean-field approximation by breaking these products. However, before under- 
taking the study of the resulting mean-field Hamiltonian, in the next section we specialize 
the above results for the CuC>2 planes in the low-temperature orthorhombic (LTO) phase of 
La 2 Cu0 4 , and discuss how to describe the structure of the corresponding classical ground 
states as commensurate spiral orders. 
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III. THE CLASSICAL GROUND STATE AS A SPIRAL ORDER 



If the rotations required at every site to bring the spin-spin interactions for all bonds in 
the lattice to the form (|3]) are compatible, then the Hamiltonian H spin = (1/2) J2(ij) is 
isomorphic to the isotropic Heisenberg model. This happens when the product of the four 
rotations around each basic plaquette in the lattice equals unity, which corresponds to the 
unfrustrated case. In this case, in correspondence with the ground-state degeneracy of @ 
there is also a large degeneracy in the ground state of H spin . ^From the relations ([|) between 
original and rotated spins we see that by judiciously chosing the direction of the Neel we 
can find ground states of H spin with no net ferromagnetic moment (S^ = — S^- parallel to d^) 
up to a maximum ferromagnetic moment of ±sin# (S' { = — perpendicular to d^). Even 
though Djf in <g) points in a special direction in spin space, because of the presence of the 
symmetric tensor T this anisotropy does not lift the degeneracy of the ground state. Q 

According to the above discussion, in order to have a net ferromagnetic moment there 
should be some degree of frustration in the system. |4j] That is, there should exist no spin 
rotations that map the entire Hamiltonian onto an isotropic one. In particular, in the LTO 
phase of La2Cu04 the crystal structure leads to the pattern of Moriya vectors shown in Fig. 
1. In this figure D 1 = Ddi and D 2 = Dd 2 , where di = (sin /3, cos /3, 0), d 2 = (cos /3, sin f3, 0), 



and D = 2Jr as in the previous section jnj. Then, the magnetic structure can be described 
in terms of two interpenetrating sublattices. Defining D 1 * 1 = (Di ±D2)/2, it has been shown 



T5| that the system develops a net ferromagnetic moment only when both D and D + are 



different from zero (f3 ^ ±7r/4). In such a case the staggered magnetization of the ground 
state is directed along D~, and the ferromagnetic moment is proportional to D + x D~ (see 
Fig. 2). 

In the two-sublattice picture the classical ground state has magnetizations = 
(S cos (ft, 0, S sin 0) and S# = (— Scoscj), 0, Ssincj)), with tan</> = rcos(/5 — 7r/4) (We are 
using the x-axis in spin space directed along D~ and the z-axis perpendicular to the Cu0 2 
plane). In order to make contact with the theory developed in it is convenient to describe 
this state as a commensurate spiral order in a distorted lattice. To this end we move the sites 
belonging to sublattices A, B by vectors 8 a = (0/vr,O) and 8b = —8a respectively (see Fig. 
3), so that the basis vectors of the new decorated square lattice are 8a and 8\ = (1, 0) + 8b 
(We took the original lattice constant a = 1 and used arrows instead of boldface to indicate 
vectors in real space). The magnetic order is then given by = (S cos Q.r, 0, S sin Q.r), 
where the magnetic wavevector Q = (it, tt) and the "Fs indicate the sites of the distorted 
lattice. This classical ground state can be described by condensing the Schwinger bosons ac- 
cording to at CT = a?a = \/~S(cosQ.r/2 + a sin Q.r/2). In the next section we will incorporate 
the quantum fluctuations on this description by considering the Schwinger-boson dynamics 
to mean-field order. 



IV. QUANTUM DESCRIPTION TO MEAN-FIELD ORDER 

Once H spin is expressed in terms of Schwinger bosons as indicated in Sec. 2 (see (|7|)), all 
the four-boson products are decoupled using as order parameters the mean values of the cor- 
relations A,B,C,D between sites in different sublattices. Then, the mean- field Hamiltonian 
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becomes: 



H MF = £[S(<f)S^- A{8)A\ - N^(\B(5)\ 2 - \A(5)\ 2 ), (8) 

i,6 5 

where = B.^+rC^ , *4 M+J - = A^-rD^ , and £(<5),.4(<f) are their corre- 

sponding mean-field values. Here N is the number of decorated unit cells and i runs on the 
A sublattice. In particular, because of the geometry of Moriya vectors in Fig. 1, the C, D 
operators take the simple form 

C ] . , = - Y ae la ^ ] a\a.^ , D. = - V e^^a^a. ,? , 

i j-i-,5 9 i. < 4<J i+o,—o ' «,i+o o ' ' Lcr i+o,<j ' 

CT Z cr 

withr ? (5 J ) = (-) J (vr/4-/3) (j = 1,4). 

We perform now a transformation to reciprocal space according to 



ft ft 

where = k+aQ/2 and the fc's are the normal modes of the decorated lattice. The resulting 
Hamiltonian has to be supplemented by adding, through Lagrange multipliers A^A^, the 
restrictions on the number of bosons on each sublattice: 



ka ka 



Notice, however, that these restrictions are only the average of the exact constraints per 
site required to have a faithful representation of the spin algebra. With this simplication 
we obtain a 8 x 8 dynamical matrix which has to be diagonalized numerically. However, 
the solution of the mean- field equations for the A, B order parameters shows that, for the 
symmetry of the spiral order, the parameters B(5) = and A^ = A#- Then, the original 
dynamical matrix decouples in two 4x4 blocks given by: 



( A - 7 (g) 1+ {q)\ 

A 7_(g) 7(g) 

-7*(<?) 7*(?) A 

V ll{q) 7*(<7) A J 



(9) 



and the corresponding one with q = k + Q/2 replaced by — q . We have defined: 
7(9) = ^E^e-^^) , 7± (g) =^r^e^ w ^ . 

Paraunitary diagonalization of (||) produces the quasiparticle dispersion relations 



<4 = A 2 - (|7| 2 + i(|7+| 2 + I7-I 2 )) ± ^j\ih + \ 2 - l7-| 2 ) Z + |7*7 + - 77 



2\ 2 _L l/v*/v «//v* 12 
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Numerical evaluation (see next section) shows that the in-plane and out-of-plane gaps are 
degenerate. For the unfrustrated case ft = n/4, so that 7± = 7, and this equation reduces 
to tJ± = A 2 — (1 +r 2 )|7| 2 . As expected, these degenerate quasiparticle bands correspond to 
an antiferromagnet with renormalized interaction J(l + r 2 ) = J / cos 2 9 = J (see (||)). 

The ground-state energy is given by E = (1/2) Z^c^A;) — (2S+ 1)N\ . The mean-field 
equations that determine the parameters A and A(5) are obtained by minimization of Eq. 
In the next section we discuss the results obtained by solving numerically these equations. 



V. RESULTS FOR THE LTO PHASE OF LA 2 CU0 4 

We have solved numerically the consistency equations by iteration, starting from the 
classical values of the order parameters and Lagrange multiplier (The classical values of the 
order parameters can be obtained by using the structure of the boson condensate given at 
the end of Sec. 3, and A c / = —E d /NS 2 = 2(1 + r 2 )). For the unfrustrated case ft = tt/4 
we obtain, as expected, two doubly-degenerate quasiparticles branches corresponding to an 
isotropic antiferromagnet with the effective coupling Jo- For ft infinitesimally smaller than 
7r/4 one of these branches develops gaps at k = 0, which correspond to the degenerate in- 
and out-of-plane gaps associated to the anisotropy. The other branch remains very much 
unaffected, i.e., almost identical to the isotropic renormalized antiferromagnetic dispersion. 
For large S the gapped branch becomes the linear spin-wave dispersion relation. P,|TTf In 
Fig. 4 we plot both branches for ft — and r = 0.4 along the typical (T,X,M) path 
there indicated. The dashed line gives the results of linear spin-wave theory for comparison. 
The fact that one of the branches remain gapless can be understood by noticing that in 
the Schwinger-boson approach the existence of a local magnetization is intimately linked to 
the boson condensation. Then, in order for the system to have a magnetized ground state, 
the chemical potential (Lagrange multiplier) must be pinned to the bottom of one of the 
original boson bands, giving a gapless quasiparticle branch. This can be confirmed by simply 
studying the isotropic antiferromagnet with a staggered magnetic field. On the other hand, 
from a physical point of view this is not a problem since the exact form of the constraint 
requires creation and destruction of both kind of quasiparticles to describe a magnon. In 
Fig. 5 we plot the g function of r for ft = [0. Notice that for the physical value 

r ~ 0.1 the result is indistinguishable from the spin wave prediction. 

The classical energy is given by E c i = —2(1 + r 2 ), and is independent of the frustration. 
In Fig. 6 we plot the ratio between the quantum and classical ground-state energies as a 
function of ft for several values of r. As can be seen, quantum fluctuations are sensitive to the 
degree of frustration, although the changes in energy are not important even for unphysically 
large values of r. In Fig. 7 we plot the behavior of the net ferromagnetic moment for the 
same values of r. As expected, |HJ it goes to zero for ft = —tt/A, where the Dzyaloshinskii 
vector D + vanishes and the ground-state configuration is completely antiferromagnetic (see 
Sec. 3). The values at ft = tt/A correspond to the ferromagnetic moments in the unfrustrated 
case (D~ = 0). However, these moments are unobservable since, as pointed out above, in 
this case the Hamiltonian becomes isomorphic to the isotropic Heisenberg model. [|J Notice 
that this happens only exactly at ft = 7r/4; for slightly different values the system develops a 
relatively large net ferromagnetic moment. Notably, as shown in Fig. 7 in the quantum case 
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the maximum ferromagnetic moment is obtained for some intermediate value of (3 which 
depends weakly on r. 



VI. INTERMEDIATE PHASE IN LA 2 -xND x CU0 4 

In the LTO phase the CuC>6 octahedra forming each CuO layer tilt in a staggered pattern 
about the (110) axis. In this section we consider the effects of DM interactions which occur 
in the presence of tilting distorsions about an arbitrary axis. Such general tilting distorsions 
may have physical relevance for La2- x Nd x Cu04. 

In order to generalize the above calculations for an arbitrary direction of the octahedra 
tilting axis we have to perform the following changes: i) In Sec. 2, the parameters t and 
6 of (Q) must be bond dependent, which implies that \Df-\ is also bond dependent; ii) 
In Sec. 3, the new pattern of Moriya vectors is given by D x = (Aj cos a, A 2 sino;,0) and 
D2 = (A2 cos a, Ai sin a, 0), with a the angle between the octahedra rotation axis and the 
[100] direction in the CuO plane, and Ai ~ (The angle a equals zero in the tetragonal 
phase and 7r/4 in the orthorhombic phase). Furthermore, the angle <fi which characterizes the 
classical ground-state spin configuration satisfies sin0 = sin if) sin 2a, where tan if) = A2/2J 



llj; iv) In Sec. 4, the parameter r = \Df-\/2J is now bond dependent. This leads to 



different A(5) for each bond direction, but 13(5) remains zero and the dispersion relations 
have the same formal expression. With these changes we obtained the following results: The 
ground-state energy and gap depend slightly on the tilting axis angle a. On the contrary, 
the net ferromagnetic moment is strongly dependent on a (see Fig. 8). It vanishes in the 
tetragonal phase (a = 0) and reaches its maximum value in the LTO phase (a = vr/4). 

In |13| it has been suggested that the inclusion of the oxygen 2p z orbitals would account 
for the weak ferromagnetism in the LTO phase of La2Cu04, with the frustration origin 
advocated in || being irrelevant. However, Entin-Wohlman et ai [11] showed that the 



effects of introducing these orbitals is too small to be responsible for the observed weak 
ferromagnetism. On the other hand, the multilevel mechanism does break the degeneracy of 
the in- and out-of-plane gaps, although the relative magnitude of this breaking is again too 
small when compared to experiments. We have also considered including the 2p z orbitals in 
our calculations. In this case the required changes are the following: i) In Sec. 2, Eq. (^) 
has to be replaced by 

Etj = JSi.Sj + Dy.(Si x Sj) + jj(l - n) [2(D y .S i )(D i7 -.S i ) - I),, 2 S„S ; 
where Q is proportional to the energy difference between the 2p z and 2p a orbitals of the 



oxygen [jlTf ; ii) In Sec. 3, the angle that describes the classical ground state has a complex 
expression in terms of the parameters of the model. We have obtained it by numerical 
minimization of the classical energy; iii) In Sec. 4, the mean-field Hamiltonian (|8|) has to 
be written as a function of the four order parameters A,B,C and D. The 8x8 dynamical 
matrix no longer decouples in two blocks, so that its paraunitary diagonalization has to be 
performed numerically. In Fig. 9 we present the new dispersion relations obtained in this 
case, for Ai = , A 2 = 0.6 , a = tc/8 , and Q = 0.4. As can be seen, the degenerate gapped 
branch splits into two branches that differ mostly near the gap wavevector. However, even 
for this large Q the ratio between the in- and out-of-plane gaps is still much smaller than 
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the experimental value. Finally, in Fig. 10 we plot the ferromagnetic moment for different 
values of Q. Remarkably, for some special value of a this moment becomes independent of 

n. 



VII. CONCLUSIONS 

We presented a Schwinger-boson approach to the Heisenberg model with DM interac- 
tion. We showed how to represent the anisotropic DM interactions in terms of Schwinger 
bosons keeping the correct symmetries present in the spin representation, which allowed us 
to perform a conserving mean-field approximation. Unlike previous studies of this model 
by linear spin-wave theory, our approach takes into account magnon-magnon interactions 
already at mean-field order and includes the effects of three-boson terms, characteristic of 
noncollinear phases, which are ignored in standard studies. Our results recover the linear 
spin-wave predictions in the semiclassical large-S 1 limit, and show a small renormalization 
in the strong quantum limit S = 1/2. 

Although the emphasis was put on the basic properties of the model Hamiltonian studied, 
for definiteness we considered a pattern of Moriya vectors corresponding to the LTO phase 
in La2Cu04. In the last section we generalized the calculations to discuss the behavior 
of relevant quantities for arbitrary directions of the octahedra tilting axis, which might 
be relevant to describe the intermediate phase in L^-^NdzCuC^. Finally, we have also 
considered the inclusion of nondegenerate 2p z oxygen orbitals, and showed that they cannot 
account for the experimentally observed difference between in- and out-of-plane gaps. 



[1] I. E. Dzyaloshinski, J. Chem. Solids 4, 241 (1958). 

[2] T. Moriya, Phys. Rev. 120, 91 {1960);ibid, Phys. Rev. Lett. 4, 228 (1960). 
[3] PW. Anderson, Phys. Rev. 115, 2 (1959). 

[4] L. Shekhtman, O. Entin-Wohlman, and Amnon Aharony, Phys. Rev. Lett. 69, 836 (1992). 
[5] S.W. Cheong, J.D. Thompson, and Z. Fisk, Phys. Rev. B 39, 4395 (1989). 
[6] J.D. Jorgensen, et al, Phys. Rev. B 38, 11337 (1988), and references therein. 
[7] D. Coffey, K.S. Bedell, and S.A. Trugman, Phys. Rev. B 42, 6509 (1990) 
[8] D. Coffey, T.M. Rice, and F.C. Zhang, Phys. Rev. B 44, 10112 (1991); ibid. 46, 5884 (E) 
(1992) 

[9] N.E. Bonesteel, Phys. Rev. B 47, 11302 (1993). 
[10] M. Crawford et al, Phys. Rev. B 44, 7749 (1991). 

[11] O. Entin-Wohlman, Amnon Aharony, and L. Shekhtman, Phys. Rev. B 50, 3068 (1994). 
[12] H.A. Ceccatto, C.J. Gazza, and A.E. Trumper, Phys. Rev. B 47, 12329 (1993). 
[13] W. Koshibae, Y. Ohta, and S. Maekawa, Phys. Rev. B 47, 3391 (1993); 48, 3580 (E) (1993). 
[14] For generality we will consider the behavior of relevant quantities for arbitrary values of (3, 

although it is in fact almost zero ||. Notice that because of the symmetries it is only necessary 

to consider — 7r/4 < (3 < 7r/4. 
[15] L. Shekhtman, Amnon Aharony, and O. Entin-Wohlman, Phys. Rev. B 47, 174 (1993). 



9 



[16] The physical values in the LTO phase of La2Cu04 are (5 ~ and r ~ 0.1 [17|. We are 

considering a larger value of r in order to emphasize the presence of a gap in this figure. 
[17] T. Thio et al, Phys. Rev. B38, 905 (1988). 
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FIG. 1. Copper sites of the CuO"2 planes in La2Cu04. We indicate the bond-depending Moriya 
vectors Dfp which alternate in sign along the x and y axis. Also shown are the corresponding 
vectors defined in the main text. 
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FIG. 3. Displacements of atoms in the A, B sublattices and the structure of the new decorated 
square lattice. 
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FIG. 4. Quasiparticle dispersion relations for r = 0.4 and (3 = 0, along the T — X — M path 
shown in the inset. The full line is the Schwinger-boson result for the degenerate gapped branch; 
the dashed line indicates the linear spin-wave prediction for comparison. The dotted line gives the 
degenerate gapless branch discussed in the main text. 
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FIG. 5. Anisotropy gap as a function of r for (3 = 0. The full line is the Schwinger-boson 
prediction; the point line gives the linear spin-wave result for comparison. 
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FIG. 6. Ratio between classical and quantum ground-state energies as a function of (3 for different 
values of r. 
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FIG. 7. Net ferromagnetic moment as a function of [3 for different values of r. 
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FIG. 8. Net ferromagnetic moment as a function of the direction of the octahedra tilting axis. 
See the main text for the definitions of the parameters involved. 
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FIG. 9. Splitting of the degenerate gapped quasiparticle branches for Q = 0.4. Upper full and 
dashed lines: Schwinger-boson results. Lower full and point lines: linear spin-wave predictions. 
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